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MODAL ANALYSIS USING A FOURIER ANALYZER, CURVE -FITTING, 

AND MODAL TUNING 

Introduction 

Since the early 1970's the dominant technique for modal testing of 
structures has been the use of single-point-excitatlon, with digital Fourier 
analysis techniques being employed for determination of frequency response 
functions (FRF's). Modal parameters (e.g. natural frequencies, damping, mode 
shapes) are derived from these FRF's by various curve-fitting techniques. 
Multlshaker sine dwell or sine sweep testing, which had predominated prior 
to the 1970*s, became less frequently employed, due to the longer test times 
and higher equipment costs Involved. A few attempts have been made to combine 
the better features of single-point FFT type testing with multlshaker testing. 
Two notable examples are References (1) and (2). 

In Reference (1) Gold and Hallauer employ a Fourier analyzer to acquire 
slngle-polnt-exdtatlon FRF's. Curve-fitting Is applied to these In order to 
determine prel Imlnary modal parameters, from which analyt’cal expressions for 
the FRF's are obtained. These analytical expressions are then employed as the 
FRF Input to the standard Asher Method of modal tuning. The Important 
characteristic of Asher's Method Is Its capability of tuning Individual modes 
In regions of high modal density. In Reference (1), the evaluation of the 
modal tuning concept using numerically simulated FRF data was very successful. 
However, evaluation In an actual modal test produced Inconclusive results 
because of the poor quality of the experimentally-acquired FRF's. 

The objective of the modal tuning procedure described In Reference (1) 
was to provide a procedure for determining accurate modal parameters (natural 
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frequencies, damping, real normal 'nodes) In situations wNere high modal density 
limits the usefulness of single-point-excitatlon techniques. Reference (2) 
also employs a tuning procedure for determining modal parameters from FRF's 
obtained by curve-fitting experimental FRF's. Whereas the original Asher 
Method^^ requires output from each shaker (Input) location and no other 
locations. Reference (2) Introduces a "minimum coincident response method" 
which permits the number of response locations to exceed the number of Input 
locations. 

In the present work a modal testing, or parameter Identification, program 
which performs the functions indicated In the flow chart of Figure 1 1s proposed, 
and portions of the proposed program are evaluated. 

Description of the Program 

The proposed modal test program alffers from single-input methods widely 
used at present In that preliminary data may be acquired using multiple Inputs, 
and modal tuning procedures may be employed to define closely-spaced-frequency 
modes more accurately or to ma'.e use of FRF's which are based on several Input 
locations. In some respects the proposed modal test program resembles earlier 
sine-sweep and s1ne-dwell testing in that broad-band FRF's are acquired using 
several Input locations, and tuning Is er<ployed to refine the modal parameter 
estimates. 

The major tasks performed In the proposed modal test program are outlined 
in the flow chart In Figure 1. They are: (1) data acquisition and FFT pro- 
cessing, (2) curve-fitting, (3) modal tuning, (4) mathematical modeling, and 
(5) computer-controlled testing. Phases (1) through (3) are described below, 
and examples are given to Illustrate and evaluate them. Phases (4) and (5) are 
the subject of further research and program development. 



Figure 1 . Major Tasks of Proposed Modal Test Program 
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Data Acquisition and FFT Processing .- This phase of the modal test program, 
which consists of acquiring FRF's based on several Input locations, resembles 
the broad>band sine sweeps commonly used In the past. However, It Is based 
on FFT processing of excitation and response records. Two forms of data 
acquisition are available: single-input and multiple-input. Figure 2 snows 
a dynamical systeiTi with m Inputs, outputs, y^(t) . The 

single-input data acquisition method Involves sequential application of a 
single shaker at each of the desired m Input locations, while the multiple- 
input method Involves operating two or more exciters simultaneously. 


X (t ) = input , y(t ) = output 

y,(t) y«(t) 

• * 0 

DYNAMICAL SYSTEM 

X,(f) “xjt) 

Figure 2. Multi pie- Input, Multiple-Output Testing 
of a Dynamical System 
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Single-Input Data Acquisitloh and FFT Processing: 

The modal tuning methods to be discussed later require FRF's which are 
based on two or more Inout locations. For example, standard Asher tuning using 
two Input locations requires that the FRF's ^21 ^22 

acquired. For minimum coincident response tuning with two Inputs and n 
outputs, the FRF's needed are H-1,2, ..., n; j«l,2). 

For single-input data acquisition, a 2-channel Fourier Analyzer can be 
used to acquire and process one FRF at a time. Preferably, all response accelera- 
tions would be acquired and recorded In a single excitation run, with post- 
processing to obtain the FRF's. 

As with other curve-fitting methods, modal tuning requires a certain 
frequency resolution In the FRF's If the tuning Is to be successful. A data 
acquisition frequency resolution on the order of one-fourth the half-power 
bandwidth Is desirable. That Is, 

Af<|(2cf^) 0) 

where r. Is a "representative (viscous) damping factor" (e.g. 4 ® 0.002) 
and Is a "representative natural frequency" (e.g. the center frequency 
of the test frequency band). Frequency resolution will be discussed In more 
detail later In conjunction with examples of Asher tuning and minimum coin- 
cident response tuning. 

Multiple-Input Data Acquisition and FFT Processing: 

Although single-input FRF's can be used In modal tuning. It Is desirable 
to acquire ♦‘he FRF's by using multiple-input?^®’®^. Consider a two-input. 
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single-output test configuration as shown In Figure 3. (The generalization 
to multiple-outputs y^(t) (1«1 ,2, . . . ,n) Is straightforward, but complicates 
the notation.) Then, 



Figure 3. Two-Input, Single-Output System 


Y(f) » Hy,(f)X^(f) + Hy2(f)X2(f) + N(f) (2) 

gives the (Fourier) transformed response, Y(f), due to the two transformed 
Inputs at Input coordinates 1 and 2. N(f) Is "noise” at the response coordinate 
which Is not linearly related to the Inputs. It can be shown^^^ that the least- 
squares estimates for H^-j and obtained by forming 



where It Is assumed that averaging Is employed and that the noise N(f) is 
not correlated with either Input. Equations (3) can be written 


®yl ‘ ”yl®n * 

®y2 “ %1®12 * V®22 


(4) 


where the 6's are the respective aver?ged auto- and cross-spectra defined 
by Equations (3) and (4). 

Assume that the ordinary coherence function between inputs x^(f) and 
Xplf) is not equal to unity. 



P 2 
1 ®22 


(5) 


That is, assume that x^ and Xg are not fuHy correlated. Then Equation (4) 
may be solved for and 
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^v1^22 


6, ^6^2 - 6,2Syi 

Hy2 ^ 


®n®22 
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Si ^12 " 

S 2 S 1 


12 


(6) 


From equations (5) and (6) It may be seen that when 6^2 * Si* “ 

A^2’ * 0» given by the slnglo-lnput expressions 





(7) 


In theory. Equations (6) may be employed to obtain FRF's when two 

Inputs are acting simultaneously, so long as the Inputs are not fully 

correlated. The above analysis can be extended to an arbitrary number 

(5 8) 

of Inputs and arbitrary number of outputs' ’ . However, If the ordinary 

coherence Is unity. Equations (6) do not hold, and a different analysis 
Is required. 

The advantages of multi-shaker testing have been discuf.sed In References 
(5) and (6). In the present situation, where modal tuning Is to be employed, 
there Is an added advantage to using simultaneous mu1t1«1nput testing to 
obtain the FRF's, since the shakers can later be employed for an actual 
tuned-dwell test based on the results of the modal tuning. 
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Curve-Fitting for r>rel1m1nary Modal Parawff tr r Estimates . - As noted by 

Gold and Hallauer*^^ useful results can be obtained by applying Asher 

tuning to analytically-synthesized FRF's rather than to the original 

(q\ 

experimentally-obtained FRF's. Preliminary studies^ ' have Indicated 
that It Is Important to Include residual ter ,,s In the curve-fitting to 
obtain parameters for use In analytically synthesizing FRF's which are to 
be used In modal tuning. 



Figure 4. Typical Freq ency Response Measurement 


Figure 4 shows a typical measured FRF and the frequency band of 
Interest, I.e. the frequency band over which It Is desired to "match" 
the measured FRF by a synthesized FRF. For general viscous damping a 
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frequency response function for Input Xj(t) and output y^(t) can be 
written^ 


H (f) . ^ fAlI_ + 1 

U Xj{f) r»l ^j2irf - s^ J^irf - s^ 




where 


V “ -®r ♦ J“dr 


th 


0 )^^ ■ damped natural frequency of the r mode 


■ decay rate of r^^ mode 


^Ijr " residue of the r^^ mode 


^IJr J'^ljr 


n ■ order of the curve fit 


Equation (8) can also be written 


'lj 


(f) ■ ; 


r»l 


r C^j^eJ^ljr C,j^e'J‘^1jr 

Jwdr - ^TTf) jo^ (a>dp + Znf) - ja 


rj 


( 9 ) 
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where C^j^, and are all real numbers. 

Although for a real system the number of modes, n , Is Infinite, 
only a limited number of modes can be onployed In the analytically-synthe- 
sized FRF. On the other hand, as shown In Figure 4, the frequency band of 
Interest may be from f^ to f|, , and In this range the number of modes 
Is finite. An expression like Equation (8) or Equation (9) can be used to 
fit the measured FRF In the frequency band of Interest, and residuals can 
be employed to approximate the contributions of modes whose frequencies lie 
below f^ or above f^^ . Then, Equation (8) can be written 


H 


u 


‘-1J * 


r«r. 



N.1r‘ \ 



( 10 ) 


where 

0 ) ■ 2irf 

r. « lower mode Index of frequency range of Interest 

O 

Tjj ■ upper mode Index of frequency range of Interest 

■ lower residual term 

Z^j ■ upper residual term (residual flexibility) 

The lower residual term can also be written^ 


L 


1j 



( 11 ) 


where 


Y 


U 


Inertia restraint 
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Thus, wUh the lower and upper residual terms approximated by use of the 
real constants (Inertia restraint) and (residual flexibility), 
the frequency response function is approximated by 



Reference(io) gives additional forms for which are equivalent to Equa- 
tion (12). For each mode Included In the summation there are four real con- 
stants. Hence, If N modes correspond to the frequency range from r. 

a 

to r^j , then Is defined analytically by 4N+2 real constants. These 

constants may be obtained by application of a curve-fit algorlthn such as 
the MDOF curve-fit algorithm In the MODAL-PLUS program. 

In the section on Examples It will be shown that obtaining analytical 
FRF’s Lv curve-fitting exper’'«ntal FRF’s prior to modal tuning serves three 
important funct1ons--it provides Initial estimates of natural frequencies, 

It permits Interpolation between experimental data points to establish more 
accurate natural frequencies, and It permits residuals to be employed so that 
"resonant nf/des" can be tuned. If desired, curve-fitting can also be used to 
generate columns of the FRF matrix corresponding to response locations where 
no physical Input was applied, This was done In Reference (2) , where a full 
FRF matrix was generated from single-point experimental FRF's. 

Mo dal Tuning for Refined Modal Parameter Estimates . - As indicated in the 
Introduction, modal tuning to be employed to refine the modal parameter 
estimates, especially where there are closely-spaced natural frequencies. 
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Figure 5 shows and KHg^) for a simple analytical 2D0F system 

with - 5.00 Hz, “ 5.01 Hz, ^ ■ 0.01, and ■ 0.01. The I(H^^) 
curve would seem to Inolcate a single mode at 5.00 Hz, while the UH^i) 

/ 5 £ \ 

curve Indicates modes at 4.98 Hz and at 5.04 Hz. Allemang' ’ * and other 
authors have proposed various means for reconciling such apparent "Incon- 
sistencies" In modal parameters based on FRF's from two or more response 
coordinates. The approach employed here Is to use modal tuning to Isolate 
the closely-spaced-frequency modes, which are frequently the source of such 
apparent Inconsistencies. 


TRANS H^^ R#i 5181 #Ai 1 EXPAND 




Figure 5. Imaginary Parts of H.j^ and H^^ for System with Modes 
at 5.00 Hz and S.Ol Hz. 
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Standard Asher Method of Modal Tuning: 

As Indicated In the Introduction, both the standard Asher Method 
Introduced In Reference 3 and a ''minimum coincident response method" 
discussed In Reference 2 will be employed In the present study. (Note: 
Reference 1 applies the term "standard Asher Method" only to the use of 
directly-measured FRF's for Asher tuning. Here we will apply the term 
more generally to tuning using the method described In Reference 3, whether 
cn directly-measured FRF's or on analytically-synthesized FRF's.) 

Application of the standard Asher Method of modal tuning begins with 
either measurement of a pxp FRF matrix, [H], over the frequency band of 
Interest, or the synthes'’^ of the [H] matrix from modal parameters obtained 
by curve-fitting the measured data. The linear relationship between the 
p output transforms and the p Input transforms Is 

lY) - [H] ih (13) 

pxl pxp pxl 

The caret indicates vectors and matrices restricted to the p input/output 

A 

locations. Next, [H] is separated into its real (coincident) and imaginary 

t\ A 

(quadrature) parts, represented by [C] and [Q] respectively, giving 

[H] “ [C] + j [Q3 (14) 

Equations (13) and (14) may be combined to give 

{?} - [C](X) + j[Q3fX) (15) 

A 

If the input Is assumed to be monophase, I.e. all components of (X) 

A 

are either In-phase or 180° out-of-phase, then {X) can be assumed to 
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be a real vector at each frequency. Then, the response has real and 
Imaginary parts 


RfY} ■ [C] (X) 


(16) 


IfY) > [Q] (X) 

(2 4) 

It has been shown' ’ that If p ■ n « the total number of degrees of 
freedom of a system, and if f^ and (X(f^)} are chosen to satisfy 


R{Y(fj,)} - [C(f^)] (X(fj,)} = {0} (17) 

nxn nxl nxl 

then each satisfying Equation (17) corresponds to an undamped natural 
fr( ,jency of the system, and each corresponding fY{f )} is an undamped 
free vibration mode shape. However, when p<n , Equation (17) becomes 

CC(f^)] {X(f^)} = fO) p<n (18) 

pxp pxl pxl 

For a nontrivial solution of this eigenproblem, it is necessary for tiie 

A 

determinant of [C] to vanish. I.e, 

det [C(f^)] • 0 (19) 

A 

Some of the roots, .will be excellent approximations to the true 
undamped natural frequencies of the system, but some will be "spurious 
frequencies." "Spurious frequencies" can be iistinguished from true 
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frequencies either by enploying several different sets of Input/response 
(41 

stations^ \ or by examining t*)e phase of the response at points other 
than the p Input points. It has also been observed^^^ that when 

A 

det [C(f)3 Is plotted versus f , steep crossings of the frequency axis 
are generally excellent approximations of true system frequencies. 
Newton's Iteration method may be employed In solving for the roots In 
Equation (19) . 




( 20 ) 


where Dj, = det[C(f|^)]. 

A 

Having Identified the approximate natural frequency f of a mode 

A <\ 

and having calculated [C(f^)] from the FRF synthesizing equation, e.g. 
Equation (12), we next determine the amplitude of the shaker forces necessary 
to tune this mode. This distribution Is calculated from Equation (18). A 

A AAA 

narrow- band sweep about f with fixed X(f) * X(f ) can be performed 

V* C 

analytically, and refined modal parameters can then be determined. 


Minimum Coincident Response Method of Modal Tuning: 

Invariably the number of response points exceeds the number of 
excitation points, even In a multi-shaker sine-sweep or sine-dwell test, 
and the modal amplitudes at many of these points may be of the same order 
of magnitude as the modal amplitudes at the excitation points. The 
minimum coincident response method Introduced In Reference (2) permits 
Information from non-input points to be accounted for In the process of 
tuning for undamped free vibration modes. 
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Let Fquatlon (13) be expanded to Include q response points, while 
the number of Input points Is p < q , and let the resulting equation be 
wri tten 


(Y) • [H]{X} - [C](X) + JCQ](X) (21) 

qxl qxp px^ qxp pxl qxp pxl 

As before, we assume that fX) Is real. Then 

RfY) • {Y)p • CC]{X) (22) 

However, since [C] Is not a square matrix. It Is not possible to determine 
unique frequencies and force appropriations In the same manner as was done 
In Equations (18) and (19). In Reference (11) Ibanez discussed this 
problem and suggested the use of a pseudo- 1nve»*se. Here, however, we will 
employ the least-squares error procedure Introduced by Ensminger and 
Turner In Reference (2). 

Let the error function be the sum of the squares of the coincident 
(real) responses, I.e. 

■ {YlJ (V)^ ■ (X)'^[C]^ [C] {iil (23) 

be minimized subject to the condition 

■ Lcj, X ■ 1 (2«) 

where 

iCjj « 1^^ row of [C] . 
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The procedure employed by Reference (2) Is to minimize the error with respect 
to the components of the Input {X} at each frequency f ; then to deter- 
mine the corresponding values of e(f) and to select those frequencies 
f which correspond to minima of c(f) as the natural frequency estimates. 
The expression for {X{f)) which minimizes e({X}, f) for a specified f 

1S<2> 

(X) ■ “t: 1 T ^ 

Equation (25) Is used to compute fX(f)) , and then Equation (23) Is 
employed to compute the resulting least-squares error. 

Program Verification 

In Figure 1 the major tasks performed In the proposed modal test 
program were outlined: (1) data acquisition and FFT processing, (2) curve- 
fitting, (3) modal tuning, (4) mathematical modeling, and (5) computer- 
controlled testing. An analytical model has been formulated and an experi- 
mental model has been fabricated for use In checking out the first three of 
these tasks. 

Experimental Model . - Since most experimental modal testing Is done on 
structures whose modal properties are unknown at the outset and for which 
there Is no validated mathematical model. It was decided to construct a 
relatively simple physical model to be used In evaluating the various phases 
of the proposed modal testing program prior to application of the techniques 
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to "real hardware." Figure 6 shows the basic elements of the lightly- 
coupled 2-beam structure constructed. The structure consists of two aluminum 
box beams 1.52m long, to which damping material (one constrained layer and 
one non-constralned layer) has been applied, and which are coupled by two tor- 
sion rods, which may be various lengths and various materials. To Increase 
the modal density, tuned "outrigger beams" are attached to the main beams. 

Both coupling torsion bars and outrigger beams are attached to the main 
beams near the node points of the fundamental free- free bending mode of the 
main beams. Figure 7 shows a single main beam with outriggers and with 
two shakers attached to the beam by long "push rods." Figure 8 shows a 
closeup of the force cell and accelerometer attachments at the end of the 
box beam. Also shown is an outrigger beam, which Is cantilevered from the 
bottom of the main beam. After several unsuccessful attempts to obtain FRF's 
with two shakers operating simultaneously, the setup was modified by Insert- 
ing a block at each end of the main beams to prevent local deformation of 
the end cross-sections. The long push rods were replaced by shorter shaker 
attachment fittings. Figure 9 shows a test configuration with two main 
beams, but with the outrigger beams removed and the shaker attachments made 
more directly. 
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Figure 7. Single Beam v<ith Vajnr 
Instrumentation 


Figure P. Driginal Force Cell and 
Accelerometer Attachments 



Figure 9. Dual-Beam with Modified Ends and Modi- 
fied Shaker Attachment 
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Instrumentation. - The Instrumentation listed In Table 1 has been employed 
for excitation of the structure, recording analog signals, performing data 
acquisition and FFT processing, and performing the computations assoclateo 
with curve-fitting, modal tuning, etc. 


Item 

Accelerometers 

Accelerometers 

Force cells 

Signal conditioners 

Analog tape recorder 

Electrodynamic shakers (50 lb.) 

Fourier analyzer 

Desktop computer 

Pen plotter 

Digital computer 


Model 

Endevco, 2265-20 
PCB, 308B 
PCB, 208-A02 
Vishay, 2120 

Hewlett-Packard (Sanborn), 3900 
MB, MB1500 

Hewlett-Packard, 5420A 
Hewlett-Packard, 9825A 
Hewlett-Packard, 7225 
Digital Equipment, PDPll/60 


Table 1. Equipment List 

For single-shaker testing, the Hewlett-Packard 5420A was used both 
to provide a random voltage signal to the shaker amplifier and to acquire 
analog force and acceleration signals and produce FRF's. For 2-shaker 
testing, the tape recorder was used to supply two pre-recorded uncorrelated 
random voltage signals to the shaker amplifiers and to record the resulting 
analog force and acceleration signals. A trigger signal was also recorded 
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on the tape so that proper phase relationships could be maintained when 
the S420A was subsequently used to produce auto>spectra and cross>spectra 
for use In computing FRF's using Equations (6). The 9825A desktop com- 
puter was used for these FRF calculations and for the modal tuning cal- 
culations, while the POPll/60 was used for curve-fitting the experimental 
FRF's. 

Examples . - Examples will be presented to Illustrate the following: calculation 
of FRF's based on dual-shaker excitation, curve-fitting of experimental FRF's, 
standard Asher tuning, and minimum phase error tuning. 

Multi-shaker Excitation: 

Tests were conducted using dual -shaker excitation of the following 
experimental models: (1) single-beam configuration as shown In Figure 7, 

(2) dual-beam configuration with outriggers, with teflon torsion bars 
connecting the two main beams, and with long shaker push rods, and (3) dual- 
beam configuration shown In Figure 9. As noted before, the analog tape 
recorder was used to record the two force channels and the two accelerometer 
channels, inese records were played back as Input to the HP 5420A, which 
was used to compute the auto-spectra and cross-spectra required to compute 
^12’ * ^12’ ^21 ^“’2 Equations 15) and (6). Data are presented 

below for the two dual-beam configurations, (2) and (3), described above. 

A dual -shaker test of configuration (2) produced the data shown 1n 
Flgures 10(a) through 10(g). Figures 10(a) and 10(b) are the two force 
auto-spectra. Figure 10(c) shows the ordinary coherence between the two 
force channels, and Figures 10(d) through 10(g) show the real and Imaginary 
parts of and , where coordinates 1 and 6 are the opposite-corner 
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shaker locations (e. 9 . see Figure 9). Although Is quite similar to 
the produce^i by a single-shaker test* It Is apparent that Is too 
*'no1sy” to be an acceptable FRF. Since the force coherence Is not equal 
to unity. Equations ( 6 ) are valid over the entire frequency range. The 
two force autospectra plots. Figures 10(a) and 10(b). Indicate that In 
the vicinity of structural resonances, the excitation levels were very low. 
However, these autospectra are quite similar to those obtained In single- 
shaker tests of the same structure. 

Several further attempts to compute FRF's based on data taken for 
configuration (2) did not produce any more acceptable FRF's than those In 
Figures 10(d) through 10(g). Several changes were then made to the experi- 
mental model, resulting In the structure shown In Figure 9. End blocks 
were Installed In the main beams to provide stiffer attachment locations 
for the force cells and accelerometers, and the aluminum shaker "push rods" 
were removed, permitting the more direct shaker attachment shown 1n Figure 9. 
Some results obtained for this test configuration are shown In Figures 11 (a) 
through 11(g). Figures 11(a) through 11(c) show the two force auto-spectra 
and the ordinary force coherence, while Figures 11 (d) through 11(g) give 
and Hg^. For comparison. Figure 12(a) Is the force auto-spectrum for 
a single-shaker test of this structure with the shaker at coordinate 1 , while 
Figures 12(b) through 12(e) are the real and Imaginary plots for the single- 
shaker FRF's and Hg^. A cursory comparison 0 ^ Figures 11 and 12 would 
seem to indicate that the single-shaker FRF's are more accurate. Since 
significant Improvement In the dual -shaker FRF's was achieved by making 
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changes In the shaker attachment fixturing, It seems likely that further 
1nH>rov(»ments in the dual -shaker FRF's could be achieved by Improvements 
of technique. 


A VB 1 M • Ml a 



Figure 10a. Force Autospectrum for Shaker No. 1 
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Figure 10b. Force Autospectrum for Shaker No. 2 
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Figures 11a, b. Force Autospectra 
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Curve-Fitting : 

The MDOF curve-fit alyorlttmi of the MODAL-PLUS program of Structural 
Dynamics Research Corpc*'ation has been used to curve-fit experimental 
FRF's. Figures 13(a) and 13(b) show the MDOF curve-fits of the dual-shaker 
FRF's of Figures 11(d) through 11(g). Table 2 lists the roots determined 
by the curve-fit algorithm. 
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Figure 13a. MDOF Curve-fit of Configuration (2). 
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Figure 13b. MDOF Curve*»f1t of Hg, . Configuration (2). 
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.819 
0.2153E-01 
104.2 
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1 .490 
•1 .223 
1 .963 
3.142 
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Phase 




"X?59 


112.5 

0.3728E-03 

0.3629E-01 

1 .474 

118.8 

0.1139E-02 

18.20 

-0.3084E-01 

119.5 

0.501 2E-02 

17.73 

-2.867 

121 .7 

0. 9861 E- 02 

1 .115 

-2.053 

129.1 

0.1089E-01 

0.2363 

-3.142 

Table 2. 

Estimated Roots 

for and Hg^ . 
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Moddl Tuning; 

As noted earlier, both the standard Asher Method of modal tunlnt^, 
which Is based on Equations (18) and (19), and a minimum coincident 
response method based on Equations (23) and (25) have been employed. They 
have been applied to the following data: (1) analytical 2D0F FRF's, 

(2) single-shaker FRF's (no curve-fitting), and (3) dual-shaker FRF's 
(no curve-fitting). Work Is currently In progress on using curve-fit 
FRF's for modal tuning. 

Figure 14 shows the simple 2D0F analytical model used to test the 
frequency separation capabilities of the standard Asher Method and the 
minimum coincident response method. Table 3 Illustrates the fact that 
the minimum coincident response method was able to separate the modes for 
the (5.00, 5.05) Hz case but not for the (5.00, 5.01) Hz case. The standard 
Asher Method, on the other hand, was able to separate the modes for the 
(5.00, 5.01 )Hz case. Table 4 shows the frequencies obtained by applying 
standard Asher tuning and minimum coincident response tuning to single- 
shaker and dual -shaker FRF's obtained for the experimental model shown in 
Figure 9. It has not yet been possible to complete the modal tuning cal- 
culations based on curve-fits of the single-shaker and dual-shaker experi- 
mental data. 
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Figure 14. 2D0F Analytical Model 



Mode 1 
S.OOHz 




Standari Asher 
Method 


5.0039 5.0156 


Minimum Coincident 
Response Method 5.0039 


Mode 1 
S.OOHz 


5.0039 5.0508 


*Not computed 


Table 3. Modal Tuning Based on 200F Analytical FRF's 


Method 


Standard 

Asher 

Method 


Minimum 

CoIncldenL 

Response 

Method 


Test Type 

Mode 1 

Single-Shaker 

119.727 

Dual -Shaker 

119.727 

Single-Sha ker 

119.922 

Dual -S ha ker 

119.727 


Mode 2 


120.117 


Table 4. Modal Tuning Based on Experimental FRF's 
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Summary and Conclusions 

A nodal testing program consisting of data acquisition and FFT 
processing, curve-fitting, modal tuning, mathematical modeling, and 
computer-controlled testing has been outlined. Both single- shaker and 
dual -shaker methods of data-acquisitlon are permitted. An experimental 
model and several analytical models have been developed for use In 
verifying steps In the modal testing program. Frequency response func- 
tions of the experimental model have been acquired using both single-shaker 
and dual -shaker methods. The particular lightly-coupled structure which 
was tested led to considerable difficulty In obtaining good FRF's using 
two shakers. Further work Is needed In order to develop multlshaker 
test procedures to the state of usefulness currently enjoyed by single- 
shaker testing. 

Computer programs have been written to Implement both standard 
Asher tuning and minimum coincident response tuning as means of refining 
modal parameter estimates through the use of actual or simulated multi- 
shaker FRF*s. Based on the application of both tuning procedures to 
analytically-formulated FRF's and to experimentally-acquired FRF's, It 
appears that the standard Asher Method may be able to separate modes In 
situations where the minimum coincident response method "breaks down." 
Further research 1s needed to determine the limitations of both of these 
tuning procedures. 

Two phases of the proposed modal test program have not yet been 
addressed. These are mathematical modeling and computer-controlled 
testing. These phases will be addressed In future research. 


39 


References 


1. R. R. Gold and W. L. Hallauer* Jr., "Modal Testing with Asher's 

Method Using a Fourier Analyzer and Curve Fitting," Proc, 25th 
International Instrumentation Symposium, The Instrument Society 
of America . May ^-10. 1979. pp. 18S-192. 

2. R. R, Ensmi/ger and M. J. Turner, "Structural Parameter Identifi- 

cation ' •om Measured Vibration Data," yAA/ASME/ASCE/AHS 20th 
Structures. Dynamics and Materials Conl^erence, St. Louis, April, 

3. Asher, G. W., "A Method of Norme’ Mode Excitation Utilizing Admittance 

Measurements," Proc. National Specialists' Meeting on Dynamics a nd 
Aeroelasticity . Ft. Worth, published by Inst, of Aeronautical 
ildence, 1966, pp. 69-76. 

4. R. R. Craig, Jr. and Y-W. T. Su, "On Multiple Shaker Resonance 

Testing," AlAA Journal . Vol . 12, no. 7, 1974, pp. 924-931. 

5. R. J. Allemang, Investigation of Some Multiple Input/Output Frequency 

Response Function Modal Analysis techniques, Fh.D. Dissertation, 
ITnlv.' of Cincinnati; ^ 

6. R. J. Allemang, et. al.. Using Dual -Input Random Excitation Modal 

Analysis. Final Report for Eqlln Structural Dynamics Laboratory, 

mr: ■ 

7. J. S. Bendat :.nd Allan G. Pierson, Engineering Applications of 

Correlation and Spectral Analysis , Chapter 8, John Wiley & Sons , 
1980. 

8. M. Richardson, "Modal Analysis Using Digital Test Systems," Scwinar 

on Understanding Digital Control and Analysis in Vibration Test 
Systems , The Shock and Vibration information Center, ilashington, 
C^TTTMay 1975. 

9. R. R. Craig, Jr,, SIR Antenna Studies , NASA Johnson Space Center, 

August, 1979 (Informal report) . 

10. D. L. Brown, et. al., "Parameter Estimation Techniques for Modal 

Analysis." SAE Transactions , v, 88, sect. 1, 1979, pp. 828-S46. 

11. P. Ibanez, "Force Appropriation by Extended Asher's Method," Paper 

760873. SAE Aerospace Engineering and Manufacturing Meeting . 1976. 



